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, Shanon's fundamental coding theorems relate classical information theory to thermodynamics. 

More recent theoretical work has been successful in relating quantum information theory to thermo- 
dynamics. For example, Schumacher proved a quantum version of Shannon's 1948 classical noiseless 
coding theorem. In this note, we extend the connection between quantum information theory and 
thermodynamics to include quantum error correction. 

There is a standard mechanism for describing errors that may occur during the transmission, 
storage, and manipulation of quantum information. One can formulate a criterion of necessary and 
\q , sufficient conditions for the errors to be detectable and correctable. We show that this criterion has 

a thermodynamical interpretation. 
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Q ■ I. INTRODUCTION 

Modern information theory was invented more than fifty years ago by Claude Shannon |33| . In his seminal paper, he 
gave a mathematical definition of information and proved his theoretical noiseless coding and noisy coding theorems. 
(— | ■ The noiseless theorem describes the physical resources necessary to store and transmit the information contained in a 
message. The noisy theorem describes the informational capacity of a noisy channel. The basic concept behind both 
theorems is that of typical messages. As one considers messages of increasing length, some messages emerge as likely 
to appear and some messages recede and become unlikely to appear. 

In statistical mechanics, one studies systems of a very large number of particles. The business of thermodynamics is 
to study the properties of the entire system that materialize as the number of particles grows. Information theory as 
conceived by Shannon can be compared to the physical science of one-dimensional statical mechanics and the concept 
■ of typical sequences is like a thermodynamic equilibrium in this statistical mechanical system. 

The success of treating information theory as a physical science is by now unquestioned and the theoretical bridge 
between information theory and thermodynamics has brought a profitable trade route between the subjects. Presently, 
quantum information theory is being hotly pursued and investigators are discovering how quantum information and 
thermodynamics are connected. One direct descendant of Shannon's noiseless coding theorem is a quantum version 
proved by Schumacher 32], which shows that entropy is important in quantum information theory. Also, connections 
between information and thermodynamics have been built using entanglement in quantum systems |l7l I28L l3l| . The 
bridge between information science and statistical mechanics lengthens to include the quantum branches of both fields 
and in this tradition we present a thermodynamic interpretation of quantum error correction. 

The concept of quantum error correcting codes has developed rapidly. First quantum error correcting codes were 
discovered in [1,13 113 ■ Error avoiding codes were discovered in [Tjj: a systematic way of building nontrivial models in 
which dynamical symmetry allows unitary evolution of a subspace [decoherence-dissipation free] while the remaining 
part of the Hilbert space gets strongly entangled with the environment. Necessary and sufficient conditions for the 
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ability of quantum error correcting codes to correct errors appeared |j, |9j, Il5l |22j, l.3(| . Let us remind the reader about 
these conditions. 

A quantum code C is a subspace of a Hilbert space W equipped with an inner product ( | ). W is sometimes called 
an encoding space. Errors are represented by a collection £ = {E a : W — > W} of linear operators. One imagines the 

diagram C c — *• ly — > PF — * C where the middle arrow represents deterioration of the medium while either during 
storage in a quantum device, or during transmission down a noisy channel. The arrow on the left represents encoding 
the information into the encoding space and the arrow on the right represents some kind of recovery procedure. The 
goal of quantum error correction is to control the code space C, the encoding procedure, and the decoding procedure, 
so that the composition of all three arrows acts as the identity on C . Here we highlight the following important result: 

Quantum Error Correction (QEC) Condition. The necessary and sufficient condition for the errors £ to be 
correctable is that (ip\ElEi,\i(}) be the same for all unit vectors ip € C and for every E a , Eb G £ . 

In the appendix, we briefly review quantum error correction and various critera equivalent to the QEC stated above. 
In particular, we show that the QEC criterion is equivalent to the widely known condition that for any orthonormal 
basis {l^j)} °f the code space C and for all E a , Eb £ £, there exists a constant c a b so that (ijjj\E^Eb\ipk) — c a b$jk- Since 
1996, several good quantum codes have been developed. Many are adapted from the classical theory of error correction 
codes, the most famous being the CSS codes H,|35| and their generalization — the stabilizer codes [E0E1E1]- See 
the textbook and the report j3|| for a summary, and also the references within. 

Let us now turn to statistical mechanics. A qubit of quantum information will be identified with a spin | and 
we consider a one dimensional system of interacting spins. The state of the system is a unit vector in a Hilbert 
space, which we will again call W. The dynamics of the system is determined by a Hermitian operator Ti : W — > W 
called the Hamiltonian. Akin to the emergence of typical sequences in Shannon's theory, or the typical subspace of 
Schumacher's theory, a special subspace of W reveals itself as the number of spins grows to infinity. Let us denote this 
subspace, which is called the subspace of thermodynamic equilibrium (see definition ^| or just thermo-equilibrium 
space, by C. 

Now, we describe some important features of the subspace C. Quantum spin chain models include a real parameter 
T called temperature. Usually, when T = 0, there is a unique vector in W, called the ground state, corresponding to 
the smallest eigenvalue of H. As a one dimensional space, the span of the ground state is too small to be used for 
quantum coding. However, for T > 0, the space of thermodynamic equilibrium has an exponentially large dimension 
(see the line preceding equation Q17[lL Our attitude is that the entire space of thcrmo-cquilibrium is a kind of 'high 
dimensional ground state'. It behaves in most ways as a single state, and any unit vector chosen from the thermo- 
cquilibrium space will serve, equally well, to represent the macroscopic physical properties of the entire equilibrium. 
Physical properties are properties such as energy, scattering matrix, and most importantly for this paper, local 
correlation functions. A correlation function is a physical quantity associated to a operator : W — > W (see equation 
(|19|l for a definition). It depends on T and is denoted by by (0)t- We present the following fundamental observation: 
For solvable models, such as XXO, XY, XYZ, Hubbard model, etc., one has [TTL ITsl Il9ll2^| : 

Thermodynamic Correlation Function (TCF) Condition. In the thermodynamic limit, the correlation func- 
tions (0)t of local operators satisfy the equation (0)t = (ip\0\tp) f or an U un it vector ip S C , the subspace of 
thermo-equilibrium. 

For any errors E a and Eb affecting an arbitrary finite number of qubits, the operator = E\Eb will be a local 
operator and the TCF conditions holds in the thermodynamic limit, thus asymptotically delivering the quantum error 
correcting conditions. The reader may wonder how one can manipulate encoded information, say for the purposes of 
computation, if no local operators act non-trivially on the code space. The TCF condition is strong, but does not 
preclude using thermodynamic systems as quantum processing devices since nonlocal operators are available. Two 
different approaches come immediately to mind. One is to use long products of Pauli operators affecting a number 
of qubits comparable to the size of the entire system. Another is to use adiabatic changes of boundary conditions to 
obtain encoded gates. We expand the adiabatic approach in the final remarks of the conclusion. 

We remark that the subspace of thermo-equilibrium is somewhat like a positive temperature, dynamic analog of 
the zero temperature ground state degeneracy that appears in top ological quantum computing — a subject that is also 
being investigated for its asymptotic ability to correct errors |12l l21j . 
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II. SHANON'S NOISELESS CODING THEOREM 



We briefly review Shannon's noiseless coding theorem. Consider a random binary variable X where the probability 
that X = is p and the probability that X = 1 is 1 — p. Consider a sequence consisting of n values of X. Then, the 
expected number of O's in the sequence is np and the number of sequences with this expected number of O's is 



(1) 



K np) (np)\(n(l — p))\ 
For large n, we have the asymptotic result 

) ~ 2 TlS(p) where S(p) = -plogp - (1 - p) log(l - p). (2) 



\npyi 



Here, and throughout this paper, log means log base 2. Note that < S(p) < 1 and S(p) — 1 only if p = |. A 
sequence 6162 • • • b n of values of X that contains the expected number of number of zeros (assuming np is an integer) 
is called a typical sequence. The probability that b\bi- ■ -b n is typical approaches 1 as n approaches infinity. So, in 
order to communicate a given sequence of length n, one needs only communicate which one of 2 nS ^ typical sequences 
is at hand. In this way n bits can be encoded in nS(p) bits. Translated into a coding theorem, one encodes blocks of 
n bits by using nS(p) bits and the probability of being able to successfully decode a block approaches one as n tends 
to infinity. 

One easily generalizes to the case that the random variable X takes values in some finite set {x}. If X takes the 
value x with probability p(x), then one defines the Shannon entropy of X to be 

S(X) = -5>(*)log(p(*)) (3) 

X 

and the same analysis and conclusions hold. The probability that a sequence X1X2 • • • x n of values of X is typical 
approaches one as n approaches infinity. 



III. SCHUMACHER'S CODING THEOREM 



A quantum bit, or a qubit, can be represented by a vector in a two dimensional complex linear space with an inner 
product. One should consider two vectors that differ by a nonzero scaler factor to represent the same qubit. An 
ensemble of n qubits is represented by an element of the n-fold tensor product of qubit spaces. Let us fix notation. 
Let W = ®" = iWj, where each Wj is a two dimensional vector space with orthonormal basis {\0)j, Each 
Wj — C\0)j © C|l)j ~ C 2 is identified with a space of qubits. The following convenient notation is suggested for a 
basis of W. Every basis vector |6x)i ® ^2)2 <S> • • • <E> \b n ) n can be referred to by the shortened name 1 6162 • ■ ■ b n ). Bit 
strings of length n index the basis vectors of W and, in this way, a general quantum state consisting of n qubits can 
be thought of as a linear combination of bit strings of length n. 

Let us define the density matrix associated to a state £ W. One has the linear functional (ip\ in the dual space 
W* whose value on \<j)) £ If is given by <fi). A matrix p can then be identified with \ip) ® via an isomorphism 
W ® W* ~ Hom(W / , W) ~ 2" x 2 n matrices. It is common to drop the tensor sign and write p = \4>)(ip\- Of course, 
the information contained in p is no different than the information contained in \tjj). If is a unit vector and {]&)} 
is an orthonormal basis of W, then p has an expansion 

p = J2p(b)\b)(b\ (4) 

b 

and the numbers p(b) are nonnegative real numbers satisfying ^2 b p{b) — 1. The numbers {p(b)} define a probability 
distribution on the collection {\b)}. One can extend the definition of p to an arbitrary \<p) £ W by p{\4>)) = (<p\p\tf>). 

Now, we can view the density matrix p as a quantum random variable which produces the state |6) with a probability 
of p(b). Or, one imagines a "quantum signal source" |32| which encodes the classical bit string b, which is the output 
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of a random variable with probability p(b), as the quantum state \b). One can then define a state I&162 ■ ■ ■ b n ) to be 
a typical state provided 6162 ■ ■ % b n is a typical sequence a la Shannon. Define the typical subspace of W to be the 
subspace T spanned by all typical states. Then the subspace T serves to compress W as the probability that a random 
vector from W will lie in the subspace T approaches 1, as n — * 00. The dimension of T is 2 nS ^ p \ where S(p) is the 
Shannon entropy defined in equation (J3J. In the present context, it is natural to note that S(p) = tr(/olog(/o)), which 
is called the von Neumann entropy of p. However, in the case that {\b}} is a collection of arbitrary, not necessarily 
orthogonal, states then it is the von Neumann entropy, not the Shannon entropy, which computes the dimension of T. 
Further refinements to the role von Neumann entropy plays in quantum compression are still being developed 



IV. QUANTUM SPIN CHAINS 



Now let us show how to formulate the quantum error correcting condition via thermodynamics. The subject of 
quantum statistical mechanics and thermodynamics is treated in several textbooks |23l |37j • There are various models 
determined by different Hamiltonians describing different interactions, and many, such as the spin chain models, have 
relevance to quantum information theory. For example, programs are underway to investigate various aspects of the 
XYZ family of quantum spin chain models with relevance to quantum computing, aspects like the implementation 
of gates pj, decoherence free subspaces 0, entanglement j^, and encoded universality Also, some two 

dimensional spin models arise in topological quantum computing [T2I l2l| . 



A. Notation 

Let a x , <r v , and a z be the Pauli matrices, which act on C 2 = C|0) © C|l) as the matrices 

It is common in quantum information theory to use the notation X, Y, and Z instead of a x , a y , and a z . Set, as 
before, W = ®? =1 W^, where Wj — C\0)j © C|l)j ~ C 2 . Sometimes the integer n is called the length of the lattice 
or the size of the model, and the reader may compare it to the length of the message in information theory. For any 
j = 1, . . . , n and superscripts a = x,y, z define an operator a" : W — > W by 

„ a , fid if i^m, 

j \w m = < a ., . (6) 
J I <t 11 j = m. 

The operator a" acts non trivially on the j-th qubit as the Pauli matrix a a and leaves all other qubits unchanged. 
One also has the matrix 

.-:=^(^-^)=(;0), (7) 

called a creation operator. By setting aj = — ^J); one can obtain every basis vector in W by acting on 

|0) := 1 00 • • • 0) by products of the aj. 

A local operator : W — > W is defined to be a linear combination of products of the operators a". If is a 
product of t such operators 

= «•••< (8) 
with distinct rrij, we say that is an operator of weight t. 
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B. The partition function and thermodynamic equilibrium 



We now wish to study the model as the spins interact. The Hilbert space of the model is W ~ (C 2 )®" and the 
interacting spins are governed by a Hamiltonian Ti : W — > W. What follows is quite general, though later we illustrate 
more details with the XXO model. So the reader may have the XXO Hamiltonian in mind: 

n 

* • ' h<T j ) (9) 

3=1 

The real parameter ft, is called the magnetic field. The XXO model with periodic boundary conditions was originally 
solved in 1961 by E. Lieb, T. Schultz, and D. Mattis [2q. The XXO model is sometimes called the "isotropic XY 
model" and is also known (in quantum information theory) as the "XY model with Zeeman splitting." 

The thermodynamic limit of the model is mathematical idealization of a very large system defined by a controlled 
limit n — ► oo. Quantities of interest are computed for finite n and then studied as n tends to infinity. Often these 
quantities are proportional to n and the proportionality factor has a finite value in the thermodynamic limit. These 
asymptotic are studied much the same way that they are in information theory when the length of the message grows 
to infinity. 

The central object of thermodynamics is the partition function Z, which is defined at a temperature T > 0, by 



Z = tr 



(2-9) . (10) 



In the thermodynamic limit, the partition function can be computed by the method of steepest descent. One has 

Z = tr2^= 2"^= J2 2" s 2~* (11) 

eigenvectors v eigenvalues E 

where the factor 2 nS is the degeneracy of the energy level E and S is entropy. Both energy and entropy increase 
linearly with n. So evaluating 



lim Z = lim y"2" 5 -T (12) 



E 

by the method of steepest descent leads to the variational equaiton 

»(*-£) -a (is) 

This brings us to the key definition: 

Definition 1. The subspace of thermodynamic equilibrium is defined to be the span of the set of eigenvectors that 
solve equation (|T3|) . 

This definition of the thermo-equilibrium subspace makes sense for all solvable models, including XXO, XXZ, XYZ, 
etc... For the model we've chosen to work with, we can be more specific. For the XXO model, equation l|12l) becomes 

Z = tr c (2^)^2- / - £iplOS ( 1+2 ^ Ei ) (14) 

where e is given by 

e(p) = -4cos(p) + 2h. (15) 

The symbol trc means the trace over the subspace C. It is defined for any operator A : W — > W by tic(A) — J>2j=i a jj 
where ajj are the diagonal entries of A when expressed as a matrix using a basis {|'0j)}| = i f° r W extending a basis 
{hfc»S=i forC. 



6 



The precise meaning of the right hand side of equation i|14|) is that linin^oo — log 2 exists and is given by 

i i r t -fen 

lim -logZ = — / cfolog(l + 2^). (16) 
n— >oo n 2ir J_n V / 

The quantity in equation is called bulk free energy. 

From equations l|34|) and (J3TJ, in section llV El we determine, asymptotically, the dimension of C. The dimension 
of C ~ 2™ s where S is the entropy: 

S = --L f dp[e(p)log0(p) + (1 - 0(p))log(l - 0(p))] , (17) 



2tt . 

where 6, called the Fermi weight, is defined by 



0(p)=(l + 2 2 * r ) \ (18) 



Physicists call e the energy of the spin wave and call p the momentum of the spin wave. 



C. Correlation functions 



Let us define correlation functions. Let be a linear combination of products of operators <r". The correlation 



function (0)t is defined by 



tr (2-T0 



(0>t = v „ ' (19) 



If has weight t, then (0)t can be related to i-point correlation function. In the thermodynamic limit, the correlation 
functions also can be calculated explicitly by the method of steepest descent, just as for the partition function 
Z = tr(2 t). Only the space of thermo equilibrium C contributes to the trace ((23, page 25). So instead of taking 
the trace over all of W, one has 



tr c 2-#0 

(0>t = ^ J - (20) 

An even stronger statement is true. In 0, 0, 0, (2^ it was shown that each term of the trace in equation H2U|) 
contributes equally, and so equation (|2Ut simplifies further: 



trc (2 



(0>t = \ ' - « for any ^ € C, (21) 

Z (ip\ip) 

which gives the thermodynamic correlation function condition stated in the introduction. 

The TCF condition is a quite general feature of quantum statistical mechanics. It holds not just for XXO, but 
also for other integrable models (XY, XXZ, XYZ, nonlinear Schrodinger, Hubbard model, etc.). We conjecture that 
the TCF condition is valid for a wide class of physically interesting models, including non-integrable models in some 
vicinity of an integrable one. 

Let us comment on the the special case T = 0. The correlation functions simplify dramatically. There is a unique 
vector \G), first identified in |26|. called the ground state. It corresponds to the lowest eigenvalue of TL. In the case 
of zero temperature 

(0>t=o = (G\0\G). (22) 



The TCF condition is precisely a generalization the equation to T > 0. 
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D. Eigenvectors of the XXO Hamiltonian 

Except for equations (|14I18|) . the discussion above applies to a quantum statistical model governed by most any 
Hamiltonian. Now, to further illustrate the thermodynamics, we work specifically with XXO. The Hamiltonian 

n 

« • n >i ■ • K) 

3=1 

can be written as H = Ho — 2hS z where 

n 1 n 

«o = - E (°Kn + and 5* = - ( 23 ) 

3=1 3=1 

Note that [Ho,S z ] — 0. The problem of finding the eigenvectors of H can be reduced to finding the common 
eigenvectors of Ho and S z . We now describe a complete set of eigenvectors of H for any positive integer n. Recall 
that aj = |(cr| - idj ). 

The eigenvectors of H are determined by a positive integer to < n and a collection of real numbers {pi\"L 1 , with 
each — tt < pi < 7r, called momenta. For short, the collection {pt}"^ may be denoted simply by {p}. Define a vector 
\{ P }) m eW by 

\{p}) m = J= J2 W {-''I' {/'!:-,, ••-,()• (24) 
The complex valued function \ is defined by 

Xm ({x}\{ P }) = ^= [] Si § n ( X& ~ *«) det ( A )' ( 25 ) 

Vto. y 1 < a<fc < m J 

where A is the to x to matrix with (j, fc) entry = exp(ixjPk) and the sign function is given by 

1 if a; > 0, 

sign(x) = { -1 if a; < 0, (26) 
if x = 0. 

The function x({ :z; }|{?'}) is called a wave function. It is symmetric in x and antisymmetric in p, so we assume that 
x\ < X2 < • • • < £ m an d that pi < P2 < ■ ■ ■ < p n - We may drop the subscript to from |{p}) m € W if the size of {p} 
does not need emphasis and just write \{p}). 

The periodic boundary conditions of the model lead to the following equation for each pj : 

exp(i Pj n) = (-l) m+1 (27) 

Note that (|27fl has n solutions and a collection {p} amounts to a choice of m of these solutions. So, for each 
to = 0, . . . , n, there are ( 71 J choices of momenta each identifying one eigenvector of H. In total, there are 2™ 
eigenvectors (|24f) as to varies from to n. A straightforward calculation of their scalar products shows that these 
eigenvectors form an orthogonal basis of W: 

(MIM) - n m and ({p}\{q}) = if {p} ± {q}. 

The eigenvalue E({p}) of |{p}) m is given by 

m 

E({p}) = J2<P^> e(p) = -4cosp + 2/i. (28) 

3=1 
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E. Thermodynamics of the XXO model 



One may think of each of the eigenvectors of TL as being obtained from the ferromagnetic state |0) by adding m 
particles (flipping m spins) via the creation operators a~ , with momenta {pj} and energies {e(pj)}. In this section we 
consider the situation as the number of particles occupying positions in the interval [— n, tt] tends to infinity (see |23jp. 
Divide the interval [— 7r, tt] into n subintervals, each of length Each subdivision point will be called a position and 
should be thought of as a possible location for a particle. They represent the n solutions of (|27|) . which the momenta 
satisfy. When a wave function has been selected, and a choice {p} of momenta has been made, one should think that 
each of the m positions appearing in {p} is occupied by a particle and the remaining n — m positions are empty (see 
figure {lj). 



FIG. 1: A picture of typical positions and momenta 
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Once a collection {p} = {pi < p 2 < 



P2 Ap 
• < p n } has been chosen, one may look at the numbers p(pj)- 

P(Pj) 



n(pj+i - Pj) 



(29) 



Now we consider the situation where the number of positions becomes large, the number of particles becomes large, 
and the density of particles remains fixed. The thermodynamic limit is defined by 



CO, 



oo, and d 



m 

— is constant. 
n 



(30) 



The reader may compare it to the double-scaling limit in matrix models. For now, consider m and n as very large, 
but still finite. As the number of particles grows, it becomes difficult to keep track of the individual momenta pj since 
there are so many of them. Instead, the growing collections {p} give way to a finite function p : [— tt, tt] — *■ R, which 
exists as a limit of expressions 1)290 . It has the interpretation that for large n 

the number of particles in the interval[p,p + Ap] ss np(p)Ap, 

provided the scale for Ap is chosen properly. One should have 

2tt 



« Ap, 



(31) 



but Ap should still be small enough to be able to approximate p by a constant on the interval [p, p + Ap] . 

One should think that p is part of a macroscopic description of the model, and that a collection {p} is part of a 
microscopic description. Be aware that each p may be the limit of many different sequences of microscopic states 
{p}, the number of which can be called the degeneracy of the macroscopic state p. For a fixed p, we can compute 
the number of microscopic states corresponding to p. First, partition the interval [—tt, tt] into subintervals of length 
Ap. It is important that the scale of Ap is such that p is well approximated by a constant on the interval [p, Ap]. 
The number of ways of selecting np(p)Ap positions from T^Ap vacancies in which to place particles is equal to the 
binomial coefficient: 



( \ 



(£Ap)! 



(32) 



\np{p)Ap) (np(p)Ap)! (n - p(p)) Ap)! ' 
Using Stirling's formula to approximate a factorial, one arrives at equation asymptotically becomes 2 nAS where 



AS = 



h log (h) ~ p{p) los(Ap)) ~ih~ p[p) ) log ih p[p) 



Ap. 



(33) 
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Now, as Ap — > 0, we find that asymptotically, p has a degeneracy of 2 nS where S, called the entropy, has the form 



S = / dp 



h log (h) - p[p) log[p{p)) ~{h~ p[p) ) log it p{p) 



(34) 



The following picture emerges. For a finite system, each collection {p} corresponds to precisely one \{p}} € W, 
which is an eigenvector of TL. In the thermodynamic limit when the number of spins becomes infinite, we replace 
the momenta {p} by the function p. But unlike the finite system, each p corresponds to many eigenvectors in the 
now infinite dimensional Hilbert space — the different microscopic states corresponding to a single macroscopic state. 
Moreover, each of these eigenvectors have the same energy density (see equation (|36|l below). The fact that every 
eigenvector corresponding to a single p has the same energy density is an instance of a crucial principle. Namely, all 
of the local observables in the thermodynamic limit depend only on the macroscopic variable p. Consider, for instance, 
the energy E({p}): 

rn rn 

E({P}) = £ e( Pj ) =n^(p j ) 2 + l~-iy (35) 
J= i J= i HPj+i Pi) 

In the thermodynamic limit, one finds that energy density is equal to 

E f 7 ^ 

lim — = / e(p)p(p)dp. (36) 

n^oo n J_ 7r 

The thermodynamic correlation function condition is a consequence of the principle stated above. The matrix elements 



(MI0IM) 



(MIM) 



(37) 



depend only on p, not on the set {p}. Therefore, in the thermodynamic limit, 



tr 2"T0 



(MI0IM) 



(MIM) 



(38) 



In order to determine which function p{p) defines the space of thermo equilibrium, we return to the variational 
equation, which arose from the steepest descent approximation. Using variational calculus, one has for XXO 



SS 



dp 



-5p(p) log(p(p)) + Sp(p) log 



2tt 



Pip) 



(39) 



and 



SE = 



dp [-6p(p)e(p)] 



(40) 



Setting 8 (S - 



nTJ 



gives the solution 



pip) = ir e (p) 

Z7T 



2vr (] 



(41) 



Therefore, the span of the eigenvectors that correspond to the function p — -^6* comprise the space C of thermody- 
namic equilibrium. Since every eigenvector in the thermo-equilibrium space corresponds to the same p, the correlation 
functions (ift\0\ijj) do not depend on which ip £ C is chosen. It can be a challenge to compute these correlation func- 
tions explicitly. In [ill Il8l . |23 | . they are computed for several models including XXO by direct calculations. 
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V. CONCLUDING REMARKS 

In statistical mechanics, the probability that an eigenstate \tp) € W will appear is given by -g2 where the 

energy E(\ip}) is the eigenvalue of the eigenvector This leads to a definition of a typical state and the subspace 
of thermo-equilibrium becomes comparable to Schumacher's typical subspace. We have shown that this subspace 
of thermo-equilibrium satisfies, asymptotically, the quantum error correction criterion for all errors of finite weight, 
thus strengthening the bond between information science and statistical mechanics. The reader may imagine that 
for finite n, an approximation of the thermo-equilibrium space may serve as quantum code. In order to develop 
this idea, one needs ways of producing gates and ma king m easurements. Gates in quantum spin chain models have 
already been implemented in several circumstances 0, l2(l l38| . Presumably, gates operating in the thermo-equilibrium 
subspace (or some finite dimensional approximation of it) can be obtained by similar means. Short products of Pauli 
matrices cannot be employed to map states in the thermo-equilibrium space into one another since such products are 
local operators and represent exactly the errors from which the thermo-equilibrium space protects against. However, 
products of Pauli matrices having a number of factors that is proportional to the length of the lattice n (the length 
of the quantum message) are candidates for gates. 

As we mention in the introduction, another promising approach for developing gates in thermodynamic codes is more 
topological. We wrote this paper assuming periodic boundary conditions. That is, the wave function Xm({£}|{p}) 
satisfies 

Xm{xi + n,x 2 , ■ ■ ■ ,x n \{p}) = Xm{xi,x 2 ■ ■ .x n \{p}). (42) 

Thus, one can imagine the lattice forming a circle. Now, by introducing a magnetic flux threading this circle, the 
boundary conditions become twisted by a real phase 4> : 

Xm(xi + n,x 2 , . . . ,x n \{p}) = e llp Xm{xi,X2 ■ . .x n \{p}). (43) 

By adiabatically changing <f> from to 2n the subspace of thermodynamic equilibrium, i.e. the codes space, will be 
mapped into itself. Even at zero temperature, this map is nontrivial — Berry's phase for such an adiabatic process was 
calculated in j24[. 

APPENDIX A: QUANTUM ERROR CORRECTION CRITERIA 

The occurrence of errors during storage or transmission of quantum data is governed by a quantum operation, also 
called a super-operator. Given a finite set of linear transformations £ = {E a : W — > W}^L X satisfying P\.Pa 
Id. one defines a super-operator Sg acting on density matrices describing states in W. The action of Sg is defined on 
a density matrix P describing states in W by 

m 

S £ (P) = J2 E ^ PE l ( A1 ) 

a=l 

Let C C W be a quantum code. One says that the code C can correct the errors £ , or that the errors £ are correctable, 
provided there exists another super-operator S-r (expressed by a collection 1Z = {Rb}t = i satisfying J2t=i ^l^b = Id) 
such that 

Sfc(Se(P)) — P f° r a U density matrices P describing states in C C W . (A2) 
We now recall (theorem 10.1 from Q) 

Theorem. Let ir : W — » C be the orthogonal projector onto the code subspace and £ = {E a : W — > WJJJLj be a 
collection of linear operators with 53^Li Pa-Pa = Id. A recovery super- operator S-r inverting Sg on density matrices 
from C exists if and only if for every E a , E S £ , there exists a constant c a b satisfying ■kE\E\ > ^ = c a i,n. 

This theorem is evidently equivalent to 
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QEC II condition. In order for the errors £ to be correctable, it is necessary and sufficient that for any orthonormal 
basis {\ipj)} l j = i of the code space C and for each E a , E <= £, there exists a constant c a b so that (ipj\E] l Eb\il>k) = c a b8jk- 

The condition QEC II as stated above seems to be familiar to those working in the held and can be found in many 
places (for example, section 3 of [22^. in chapter 7 (pages 9 and 86) of (2^, and section 6.4 of Q.) In the introduction, 
we stated the 

QEC condition. The necessary and sufficient condition for the errors £ to be correctable is that (ip\E^Eb\ip) be the 
same for all unit vectors ip € C and for every E a , Eb € £. 

We now prove that the conditions QEC and QEC II are equivalent. 

Proof. Suppose that {tp\E a Eb\vb) = c a b for every unit vector \ip) £ C. Note that the matrix A = [c a b] is Hermitian, 
hence there exists a unitary U = [u rs ] such that UAW = D = [d a b] is a real diagonal matrix. Now let us dchne 
F r = X)a=i u ra,E a - Note that for any unit vector ip £ C , we have 

&\F}F s \rJj) =d rr o rs . (A3) 

The operator F^F r is Hermitian so we can find an orthonormal basis {|</^)}j = i for C of eigenvectors of F^F r . By 
substituting \<pj) for in equation IA3I one finds that the eigenvalue corresponding to \4>j) is d rr for every j. Thus, 
F^F r — d rr I . It follows that for any orthonormal basis {lipj}}^^ of C 

(iP j \F}F s \i> k }=d rr S rs 8 jk . (A4) 

By changing back from the {F r } to the {E a } by using E a — J2™=i u r aF r , we find that 

(^•l^blVfc) =c ab 6 jk . (A5) 

This proves that QEC QECII. 

Now, suppose that {|*0j)}j=i be an orthonormal basis for C and that (iftj \E\Eb\4jk) — c a bSj k - For any unit vector 

ijj e C we have ip = Ylj=i a j*Pj f° r some a 3 with Ylj=i °^j a j = 1- We compute 

1 1 

(mim) = y^zm*s\Kei,W) (A6) 
3=1 /c=i 
1 

3=1 

= c ab . (A8) 
Thus QECII QEC. In such a way we proved the equivalence of QEC II and QEC. □ 
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